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Abstract: Cross sections for the production of pairs of photons plus two additional jets produced from double parton 
scattering in high-energy proton-proton collisions at the LHC are calculated for the first time. The estimates are 
based on the theoretical perturbative QCD predictions for the productions of 77 at next-to-next-to-leading-order, jet 
-h jet and 7 -I- jet at next-to-leading-order, for their corresponding single-scattering cross sections. The cross sections 
and expected event rates for 77 -|- 2 jets from double parton scattering, after typical acceptance and selections, are 
given for proton-proton collisions with the collision energy ^/s = 13 TeV and an integrated luminosity of 100 fb~^ 
planned for the following two years, and also ^/s = 14 TeV with 3000 fb~^ of integrated luminosity as LHC designed. 
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1 Introduction 

In proton-proton (pp) collisions with higher energies 
at the Large Hadron Collider (LHC), particle production 
is dominated by multiple interactions of their constituent 
partons, with most particles from the hardest proton- 
proton scattering and the radiation and fragmentation 
of secondary partonic actions. The higher centre-of-mass 
energy leads to enhanced parton densities which cause a 
sizable probability for two or more parton-parton scatter¬ 
ings within the same pp interaction [H, At LHC, vari¬ 
ous measurements of the differential distributions in W -I- 
jets [1,13 and J/tp + W i show that the excesses above 
the expectations from single parton scattering (SPS) are 
consistent with double parton scattering (DPS). Vari¬ 
ous measurements in other pp and pp collisions at ^/s 
— 63 GeV Hi, 630 GeV Q, and 1.8 TeV H are consis¬ 
tent with DPS contributions to multijet final states, as 
well as to 7 -I- 3-jet events at ^/s = 1.8 TeV H and 1.96 
TeV 0 - The measurements of DPS processes can pro¬ 
vide valuable information on the transverse distribution 


of partons in the proton n| and on the multi-parton cor¬ 
relations in the hadronic wave function [ll|. DPS also 
constitutes the background for new physics searches at 
the LHC 13- 13 . Additional searches for DPS have been 
proposed via double Drell-Yan, four jets and same-sign 
WW production 16l-[T8|. 


double parton scattering (DPS) in high-energy proton- 
proton collisions at the LHC are calculated for the first 
time. 77 final states have played a crucial role in the 
recent discovery of a new boson at the LHC [3 113 and 
are also important in many New Physics searches [El- 
24| , in particular the search for extra spatial dimensions 


In this paper, the cross section for the production of 
pairs of photons plus two additional jets produced from 


or cascade decays of heavy new particles. In particular, 
diphotons in combination with jets and missing energy 
occur in gauge mediated SUSY scenarios. 77 or plus 
two additional jets offers also an important test of both 
perturbative and n on-p erturbative quantum chromody¬ 
namics (QCD) [25l427| . 77 - 1-2 jets is also the main 
irreducible background for other physics analyses with 
77 and jets in the final state at the LHC, such as Higgs 
produced in vector boson fusion. For the production 
of 77 -I- 2 jets, a sizeable contribution from DPS with 
77 produce in one scattering while the second scattering 
yielding two jets can be expected. 

The structure of this paper is organized as follows. In 
section 2, a generic way of the DPS cross section as the 
product of the SPS cross sections and its parameter are 
briefly introduced. The details of the cross section of 77 
-I- 2 jets calculation and the cross section of different SPS 
processes estimated from higher order theoretical predic¬ 
tions are described in section 3. The results including the 
cross section of 77 -|- 2 jets with y/s = 13 TeV and 14 TeV 
at LHC and the expected event rates, with typical selec¬ 
tions,a re summarized in section 4. The summary and 
outlook are given in section 5. 
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2 Generic formule of DPS 


For a composite system {A-\-B) in hadronic collisions, 
its production cross section from DPS, can be 

written model-independently as the product of the cross 
sections of A and B originated from single parton scatter¬ 
ing, and normalized by an effective cross 

section a^ff [28 1 


77 - 1-2 jets from DPS with ^/s = 13 TeV and 14 TeV at 
LHC. A number 5 mb was assigned as its uncertainty to 
estimate its effects on the final results. The uncertainty 
on (Te// is the dominant uncertainty for the calculation 
of the production cross section of 77 -I- 2 jets from DPS 
in the following sections 

3 (T^p^^^+ 2 jets calculation 


OPS _ 

^pp^AB ~ 


SPS „ SPS 

m Ppp^B 
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According to the descripition in above section, the 
production cross section of 77 -I- 2 jets from DPS in pp 
collisions can be written as 


where m ia a symmetry factor accounting for distinguish¬ 
able (to= 2 ) and indistinguishable (m=l) final-states. 

The effective cross section is a measure of the 
transverse distribution of partons inside the colliding 
hadrons and their overlap in a collision. It is indepen¬ 
dent of the process and of the phase-space under consid¬ 
eration. A number of measurements of Ue// have been 
performed in pp and collisions at ^/s = 63 GeV 
630 GeV Q, 1.8 TeV [a. 1^, [ 2 ^ . 1.96 TeV [lC)| and also 
7 TeV at LHC ii- The measured values range from 
5mb at the lowest energy to about 20mb from CMS at 
7 TeV. Fig. □ shows a comparison of the effective cross 
section / measured by different experiments using dif¬ 
ferent processes at various centre-of-mass energies. 
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Fig. 1. Ue// measured by different experiments 
using different processes i, 0^13, Ei]. The 
’’Corrected CDF” data point indicate the Ue// 
value corrected for the exclusive event selection 


■CMS (W + 2jets) 

■ ATLAS (W + 2 jets) 

- CDF (4 jets) 

- CDF (y + 3 jets) 

- Corrected CDF {y + 3 jets) 

- DO {y + 3 jets) 

■ UA2 (4 jets - lower limit) 

■ AFS (4 jets - no errors given) 


The measured values of a^/f from TeV experiments 
at Tevetron (CDF and DO) and LHC (ATLAS and CMS) 
are consistent with each other within their uncertainties. 
In the following calculations, a numerical value Ue// ~ 15 
mb was used to estimate the production cross section of 


DPS 

^ pp — 
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pp—^ pp—¥2jets 
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pp—¥^-\-jet 


X cr; 


pp — 
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77 production has been calculated at next-to-leading- 
order (NLO) some time ago [s^, supplemented also by 
gluon initiated subprocesses beyond the leading order 
311 and soft gluon resummation [^, [^. Recently, 
next-to-next-to-leading-order (NNLO) corrections to di¬ 
rect diphoton production also have become available 341. 
The measurements from LHC (25l - [27| show that the 
NNLO can give much better agreement with measured 
data than the lower order predictions. For the integrated 
cross section, the predicted values by NNLO are almost 
exactly the same 27j or consistent within the uncertain¬ 
ties 25, 2^ with the measured ones. So for the produc¬ 


tion cross sections of 77 final-state from SPS, 
with different -/s will be obtained from the NNLO cal¬ 
culation with the package 27 NNLO. 

For the dijet cross section, the measured data at LHC 
[3l,[3^ can be well described by NLO perturbative QCD 
(pQCD) calculations from NLOJet-|—I- program [s^ cor¬ 
rected to account for non-perturbative and electroweak 
effects. From 3^, the non-perturbative correction is 


within 3% for jet reconstructed with the anti-fct cluster¬ 
ing algorithm and distance parameter or cone size 
R=0.4. The corrections for the electroweak effect can be 
negligible if the dijet mass less than about 1 TeV. In this 
analysis, the NLO calculations of cF^pS 2 jets performed 
using NLOJet-|--|- (version 4.1.3) within the framework 
of the fastNLO package (version 2.3.1) [s^. 

NLO pQCD prediction from the program JETPHOX 
(version 1.3.1) |40| is used for the calculation of j + jet 
cross sections from SPS in this paper. This program 
includes a full NLO QCD calculation of both the direct- 
photon and fragmentation contributions to the cross sec¬ 
tion. The number of flavours was set to five. Compared 
with the measurements of 7 -b j et cross sections at the 
LHC, the predictions from JETPHOX multiplied by a 
factor close to unity for the corrections of hadronisation 
and underlying-event effects give a good description of 
the Ay and measured cross sections [mEi. 
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Different PDF sets are used for the calculations of 
these three SPS processes. MSTW2008NNLO [i^ is 
used for calculations with 27 NNLO. CTIONLO 


.SPS 


^2jets 


with NLOJet++ in 


44j is used for both 

fastNLO package and with JETPHOX. 


pp—^7+jet 

The calculations of are performed with the 

factorization and renormalization scales equal to the in¬ 
variant mass of two photons, fip = fin = . The scale 

uncertainty and PDF uncertainty are also considered. 
A simplified and less computationally intensive estimate 
of the renormalization {fj-p) and factorization (/Tf) scale 
uncertainties is performed by varying these scales simul¬ 
taneously by a factor of two up and down around 
^J-F = ^J'R = and 41 eigenvec¬ 

tor sets of MSTW2008NNLO are used to build the PDF 
uncertainty envelope. 

Calculations of (Jpp^ 2 jets derived using NLO- 
Jet-|--I- within the framework of the fastNLO package 
at a factorization and renormalization scale equal to the 
average transverse momentum {p^") of the two jets {pp 
= Mh = The uncertainty due to the choice of 

factorization and renormalization scales is estimated as 
the maximum deviation at the six points (pip/p, Pr/ p)= 
(0.5, 0.5), (2, 2), (1, 0.5), (1, 2), (0.5, 1), (2, 1) with p = 
Py”®. 52 eigenvector sets of CTIONLO are used to build 
the PDF uncertainty envelope. 

For NLO calculations of using JETPHOX, 

the renormalization, factorization and fragmentation 
(pf) scales are chosen to be photon’s transverse momen¬ 
tum, Pp = Pp = pf = Ep. Same as above, 52 eigenvector 
sets of CTIONLO are used to build the PDF uncertainty 
envelope. 

Above calculations were performed with the strong 
coupling constant at two-loop order with a^imz) = 0.118 
in CTIONLO and 0.117 in MSTW2008NNLO. The un¬ 
certainty on as{mz) was not considered in this study. 
The uncertainty from scales, pdf and ajmz) is around 
10%, 5% and 1% respectively 25-13, EH, [Hi, [U, l42j| . 
Compared to the larger uncertainty on the a^ff with 
more than 30% used in this study as told in the end of 
section 2 , the effect on the final results from the uncer¬ 
tainty on as(mz) can be negligible. 


4 Results of expected 

event rates at LHC 


In this paper, several sets of typical selections at LHC 
were tried to calculate the production cross section of 77 
-I- 2 jets from DPS in pp collisions, cFpp'/j.^.^p 2 jeta- Due to 
the high level trigger requirements for 77 events at LHC 
for the higher energy and higher luminosity collisions, 
five sets of requirements on the photons’ transverse mo¬ 
mentum were considered, (EJ% E/?) > (30, 20) GeV, 
(30, 30) GeV, (40, 20) GeV, (40, 30) GeV and (40, 


40) GeV with 71 representing the maximum Ep photon 
and 72 the minimum Ep one of two photons. So for sin¬ 
gle photon requirement in the 'j+jet, 3 cases with Ep > 
(20, 30, 40) GeV were considered. The photon should be 
also constrained in the pseudorapidity region \ri\ <2.5. 
An isolation requirement is applied on the photon to 
fulfill the isolation requirement from experimental mea¬ 
surements 23-27, 41, 421. The standard isolation, the 
Ep sum of partons in a cone of size Ai?=0.4 around the 
photon required to less than 5 GeV, is applied in JET¬ 
PHOX for the calculation of Tor 27 NNLO, 

the smooth Frixione isolation [43] on the photons is ap¬ 
plied 


^ ^ ^ Vl-cos(Ai?o)^ 


(3) 


where Ep°(AR) is the Ep sum of partons in a cone of 
size AR, ARq = 0.4, e = 5GeV, and n = 0.1. This 
criterion is found to have the same efficiency as the stan¬ 
dard isolation used for the other generators within a few 
percent 23, 23]. Additional the angular separation be¬ 
tween two photons is required to be at least larger than 
0.4 {ARy.f >0.4) to ensures one photon will not enter 
the isolation cone of the other photon, which is similar 
to the requirement applied in the data analyses at LHG 
experiments ATLAS and CMS 23- 27|. 

In this study, jet is reconstructed with the anti-fct 
clustering algorithm and cone size R=0.5. Jets are in 
the acceptance region with \r]^‘'*' \ <4.5. Two tries on jet 
p^P^ were performed, pt/:^ > 20 or 25 GeV. For the dijet 
events, two jets should be separated by requiring their 
angular distance ARjj greater than 1.0 to avoid the over¬ 
lapping of two jet cones. For the 7 -I-jet production, the 
angular distance between 7 and jet should be greater 
than 0.5 (AR^j > 0.5) to ensure that the partons belong 
to the jet will not enter to the isolation cone of 7 . 

Fig. [2] shows the cross sections of computed 

from the 27 NNLO at y/s = 13 TeV and 14 TeV with 
scales and pdf uncertainties considered, for different sets 
of Ep requirements on diphotons. The scale uncertainty 
is around 10% and the pdf uncertainty is about 4%. The 
selection sets in x-axis are the 5 sets of requirements on 
(El\ EJ?), number 1 for (E?% E/?)> (30, 20) GeV, 2 
for (A?, E]?)> (30, 30) GeV, 3 for (E/j, E/?)> (40, 
20) GeV, 4 for {E/,\ E/?)> (40, 30) GeV and 5 for {E/,\ 
Ep^)> (40, 40) GeV. The detailed values can also be 
found in Table [TJ For the central values, the cross section 
with y/s = 14 TeV is about 9% higher than that with y/s 
= 13 TeV with the same selection requirements, which 
is within the scale and pdf uncertainties. 
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Fig. 2. Predicted cross sections of 77 by 2'y NNLO 
at ^/s = 13 TeV and 14 TeV with selection set 1 
for (FIJI, £;72)> ( 3 Q^ 20) GeV, 2 for (EJJ F;J=)> 
(30, 30) GeV, 3 for (SJi, E:p)> (40, 20) GeV, 
4 for (FIJI, ( 4 Q^ 3 Q) Qgy and 5 for (FIJJ 

F1J^)> (40, 40) GeV. Scale and pdf uncertainties 
are included. 


Table 1. Gross sections in unit of pb of 77 pre¬ 
dicted by 27 NNLO at ^/s = 13 TeV and 14 TeV. 
The uncertainties include the scale and pdf un¬ 
certainties. 


(Eji, EJ?)> 


13 TeV 


14 TeV 

(30,20) GeV 

89.3 

± 

9.7 

96.9 

± 

10.5 

(30,30) GeV 

44.9 

± 

4.9 

48.7 

± 

5.3 

(40,20) GeV 

48.0 

± 

5.2 

52.1 

± 

5.7 

(40,30) GeV 

31.2 

± 

3.4 

33.7 

± 

3.7 

(40,40) GeV 

18.0 

± 

2.0 

19.6 

± 

2.1 


Fig. |3] shows the differential cross sections, as a func¬ 
tion of the Pt of leading jet with both jet pj"* > 20 GeV 
and \r]\ < 4.5, of <J^p^ 2 jet 3 computed from NLOJet-|—I- 
within the framework of the fastNLO package at y/s 
= 13 TeV and 14 TeV with scales and pdf uncertain¬ 
ties plotted in the same figure. The bottom two plots 
show the relative uncertainties including the scale uncer¬ 
tainty and scale©pdf uncertainty combined in quadra¬ 
ture. The contribution of pdf uncertainty is tiny. The in¬ 
tegrated cross section are 117.6jJo(scale)^^^(pdf) (/rb) 
and 122 . 315 ®(scale) - 1 . 4 (pdf) (^tb) for y/s = 13 TeV 
and 14 TeV with both jet > 20 GeV and < 
4.5, 52.2lj4(scale)lo;g(pdf) (/rb) and 56. 2 I 2 . 3 (scale) 
lg;g(pdf) (fih) for y/s = 13 TeV and 14 TeV with both jet 
> 25 GeV and < 4.5. When pt requirements 
on both jet increase 5 GeV from 20 GeV to 25 GeV, the 
cross section are reduced almost by a factor of 2 . 



Fig. 3. Differential cross sections of jet+jet com¬ 
puted with NLOJet-b- 1 - within the framework of 
the fastNLO package, after the requirements on 
both jets with pj®* > 20 GeV and < 4.5. 
The solid black circles are the results for y/^ — 
13 TeV while the red diamond for y/s = 14 TeV. 
Bottom two plots show the relative uncertainties 
including the scale uncertainty and scale©pdf un¬ 
certainty with y/s = 13 TeV in the middle plot 
and y/s = 14 TeV in the bottom plot. 


Combined the photon FiJ requirements for the 77 
productions and the jet * requirements for the jet+jet 
productions, cross section of with six sets of 

selections on the transverse momentums of photon and 
jet with (F;J, pj"*) > (20, 20) GeV, (30, 20) GeV, (40, 
20) GeV, (20, 25) GeV, (30, 25) GeV and (40, 25) GeV 
were calculated. Fig. |4] shows the differential cross sec¬ 
tions of as a function of photon FiJ with (FiJ, 

Ft*) > (40, 20) GeV, < 2.5, \ < 4.5 and separa¬ 

tion AR^j > 0.5, computed from JETPHOX with y/s = 
14 TeV. The contributions from the direct photon pro¬ 
duction and photon from fragmentation are shown in the 
same plot. The scales and pdf uncertainties are also plot¬ 
ted in the same figure. The integrated cross sections are 
listed in Table [2] for different sets of selections and col¬ 
lision energies. The scale and pdf uncertainties are also 
listed in this table, with about 10 % uncertainty from 
scales and around 4% from pdf. 
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(GeV) 

Fig. 4. Differential cross sections of 7 + jet as 
a function of the photon E^, computed with 
JETPHOX and the selections (EJ, > (40, 
20) GeV, |ry'’'| < 2.5, < 4.5 and separation 

t^Rjj > 0.5 at ^/s = 14 TeV. The solid circle 
are the total contributions while the blue triangles 
represent the direct contribution and the red tri¬ 
angles are the contributions of photon from frag¬ 
mentation. Scales and pdf nncertainties are also 
shown in this plot. 

Table 2. Cross sections in unit of 10® pb of j+jet 
predicted by JETPHOX at y/s — 13 TeV and 
14 TeV. The total uncertainties including scale 
uncertainty and pdf uncertainty are also list in 
this table. 


{E/, 

y/s = 13 TeV 

yT = 14 TeV 

(20.20) GeV 

(30.20) GeV 

(40.20) GeV 

(20.25) GeV 

(30.25) GeV 

(40.25) GeV 

20 . 2/11 

85.0+|-^ 

97.211?:^ 

52.7/1/ 

22 . 0 / 1 / 

91.7+l°f 

45.011:1 

21.411:1 


According to Eql2]and the above cross sections of the 
SPS processes, the cross sections for the production of 
pairs of photons plus two additional jets produced from 
double parton scattering (DPS) in high-energy proton- 
proton collisions at the LHC are calculated for the first 
time. The results are summarized in Table [3l Two jets in 
the same pp — >■ 'y"f+2jets event from DPS have the same 
Pt cut thresholds, bothp^*"* > 20 GeV or 25 GeV simulta¬ 
neously. The calculated cross section can be around 0.1 
pb to «1 pb with the selections considered in this paper. 
The uncertainty on the cross section ia around 50%, with 
the dominant contribution from the uncertainty of (Te//. 

Table 3. Cross sections in unit of pb of 
iypp^-y-y+ 2 jets calculated for y/s = 13 TeV and 


14 TeV with the selections described in the pa¬ 
per. The total uncertainties including scale un¬ 
certainty, pdf uncertainty and also the Ue// un¬ 
certainty are also list in this table. 


, both 

y/s = 13 TeV 

^/s = 14 TeV 

(30.20.20) GeV 

(30.20.25) GeV 

(30.30.20) GeV 

(30.30.25) GeV 

(40.20.20) GeV 

(40.20.25) GeV 

(40.30.20) GeV 

(40.30.25) GeV 

(40.40.20) GeV 

(40.40.25) GeV 

0.84611:111 

0.428l°:ll® 

0.43111111 

0.428l°'il® 

u.‘io< _o .222 
n ooQ"i"0-113 

n 977+0137 
<- 0.141 

0 136+0-068 
U.iOD_Q Qgy 

n 1 fv4+0-076 
U.iD4t_Q Qyg 

n 076+0-038 

U.U/0_Q Q3y 

u.you_Q 47g 

U.OUU_Q 250 

0.489lO-2« 

0.2501°:121 

u.4yo_Q 247 

0 261 +0-130 
U.ZDi_o.i29 

0.313l°:l|| 

0 159+0-078 

u.±oy_Q 078 

0.176tg;X 

0.089l°:°14 


With an integrated luminosity of 100 fb~^ at y/s = 

13 TeV accumulated in the following years, about 85k 
pp — 77 -|- 2 jets events from DPS can be obtained with the 
loosest selections, diphoton (AJ% E'^)> (30, 20) GeV 
and both jets p^^ > 20 GeV. These events can be trig¬ 
gered by the diphoton paths proposed at the LHC for 
y/s = 13 TeV. When the integrated luminosity increas¬ 
ing at yfs = 14 TeV, tighter Et thresholds on diphoton 
for the trigger will be used. With the tighter selections, 
diphoton {E'jp-, E'^)> (40, 30) GeV and both jets > 
20 GeV, about 940kpp—7> 77 - 1 - 27 ets events from DPS can 
be obtained with an integrated luminosity of 3000 fb~^. 
Even with the tightest selections studied in this paper, 
diphoton (EJ\ E^^)> (40, 40) GeV and both jets > 
25 GeV, we can also get about 260k pp ^ 77 + 2,jets 
events from DPS with 3000 fb~^ as designed by LHC. 

5 Summary and outlook 

In this paper, the cross sections for the production 
of pairs of photons plus two additional jets produced 
from double parton scattering in high-energy proton- 
proton collisions at the LHC with y/s = 13 TeV and 

14 TeV (LHC Run2) are calculated for the first time. 
With the generic formula, the cross sections have been 
computed based on the theoretical perturbative QCD 
predictions for the productions of 77 at next-to-next-to- 
leading-order , jet -I- jet and 7 -I- jet at next-to-leading- 
order, with their corresponding single-scattering cross 
sections. From the LHC measurements with the collision 
data obtained in years 2011 and 2012 (LHC Runl), these 
theoretical predictions for these three SPS processes can 
give the best agreements with the measured data. With 
the typical acceptance and selections used at LHC, the 
cross sections cr//l/yj+ 2 jets can be estimated to be around 
0.1 pb to 1 pb with the collision energy y/s = 13 TeV 
or 14 TeV. The expected event rates for 77 -|- 2jets from 
DPS, with some sets of selections, are given for proton- 
proton collisions with the collision energy y/s = 13 TeV 
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and an integrated luminosity of 100 fb~^ planned for the 
following two years, and also y/s = 14 TeV with 3000 
fb~^ of integrated luminosity as LHC designed. The un¬ 
certainties on the cross section and the events rates are 
mainly dominated by the a^ff uncertainty. The scale 
and pdf uncertainties for the productions of these three 
SPS processes are also considered. 

With the incoming LHC Run2 data, there are enough 
pp —>■ 77-|-2jets events from DPS for investigations. It 
needs further studies on the variables, such as the angles 


between two photons and two jets, to be chosen for the 
discrimination of pp— -\-2jets events from DPS and 
pp —>■ + 2jets events from SPS when performing the 
data analysis. Also the contributions from the DPS to 
the whole pp —>■ ^^+2jets event rates on the distributions 
of some typical variables are need detailed investigations 
in the LHC Run2 data analysis. 

The authors would like to thank Dr. Hua-Sheng Shao 
from CERN for helpful discussions. 
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